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COLLISIONAL DARK MATTER AND THE STRUCTURE OF DARK HALOS 
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ABSTRACT 

We study how the internal structure of dark halos is affected if Cold Dark Matter particles are assumed 
to have a large cross-section for elastic collisions. We identify a cluster halo in a large cosmological N- 
body simulation and resimulate its formation with progressively increasing resolution. We compare 
the structure found in the two cases where dark matter is treated as collisionless or as a fluid. For the 
collisionless case the overall ellipticity of the cluster, the central density cusp and the amount of surviving 
substructure are all similar to those found in earlier high resolution simulations. Collisional dark matter 
results in a cluster which is more nearly spherical at all radii, has a steeper central density cusp, and has 
less, but still substantial surviving substructure. As in the colisionless case, these results for a "fluid" 
cluster halo are expected to carry over approximately to smaller mass systems. The observed rotation 
curves of dwarf galaxies then argue that self-interacting dark matter can only be viable if intermediate 
cross-sections produce structure which does not lie between the extremes we have simulated. 

Subject headings: dark matter - galaxy: formation - methods: numerical 



1. INTRODUCTION 

Cold dark matter scenarios within the standard infla- 
tionary universe have proved remarkably successful in fit- 
ting a wide range of observations. While structure on 
large scales is well reproduced by the models, the situa- 
tion is more controversial in the highly nonlinear regime. 
Navarro, Frenk & White (1995, 1996, 1997; NFW) claimed 
that the density profiles of near-equilibrium dark halos can 
be approximated by a "universal" form with singular be- 
haviour at small radii. Higher resolution studies have con- 
firmed this result, finding even more concentrated dark 
halos than the original NFW work and showing, in addi- 
tion, that CDM halos are predicted to have a very rich 
substructure with of order 10% of their mass contained in 
a host of small subhalos (Frenk et al 1999, Moore et al 
1999a, 1999b, Ghigna et al 1999, Klypin et al 1999, Got- 
tloeber et al 1999, White & Springel 1999). Except for 
a weak anticorrelation of concentration with mass, small 
and large mass halos are found to have similar structure. 
Many of these studies note that the predicted concentra- 
tions appear inconsistent with published data on the rota- 
tion curves of dwarf galaxies, and that the amount of sub- 
structure exceeds that seen in the halo of the Milky Way 
(see also Moore 1994; Flores and Primack 1994; Kravtsov 
et al 1998; Navarro 1998). 

It is unclear whether these discrepancies reflect a fun- 
damental problem with the Cold Dark Matter picture, or 
are caused by overly naive interpretation of the observa- 
tions of the galaxy formation process (see Eke, Navarro & 
Frenk 1998; Navarro & Steinmetz 1999; van den Bosch 
1999). On the assumption that an explanation should 
be sought in fundamental physics, Spergel & Steinhardt 



(1999) have argued that a large cross-section for elastic 
collisions between CDM particles may reconcile data and 
theory. They suggest a number of modifications of stan- 
dard particle physics models which could give rise to such 
self-interacting dark matter, and claim that cross-sections 
which lead to a transition between collisional and collision- 
less behaviour at radii of order 10 - 100 kpc in galaxy ha- 
los are preferred on astrophysical grounds. Ostriker (1999) 
argues that the massive black holes observed at the cen- 
tres of many galactic spheroids may arise from the accre- 
tion of such collisional dark matter onto stellar mass seeds. 
Miralda-Escude (2000) argues that such dark matter will 
produce galaxy clusters which are rounder than observed 
and so can be excluded. 

At early times the CDM distribution is indeed cold, so 
the evolution of structure is independent of the collision 
cross-section of the CDM particles. At late times, however, 
a large cross-section leads to a small mean free path and 
so to fluid behaviour in collapsed regions. In this Letter 
we explore how the structure of nonlinear objects ("dark 
halos") is affected by this change. We simulate the forma- 
tion of a massive halo from CDM initial conditions in two 
limits: purely collisionless dark matter and "fluid" dark 
matter. We do not try to simulate the the more complex 
intermediate case in which the mean free path is large in 
the outer regions of halos but small in their cores. If this 
intermediate case (which is the one favoured by Spergel & 
Steinhardt (1999) and by Ostriker (1999)) produces non- 
linear structure intermediate between the two extremes we 
do treat, then our results show that collisional CDM would 
give poorer fits to the rotation curves of dwarf galaxies 
than standard collisionless CDM. Further work is needed 
to see if this is indeed the case. 
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2. THE N-BODY/SPH SIMULATION 

Our simulations use the parallel tree code GADGET 
developed by Springel (1999, see also Springel, Yoshida & 
White 2000b). Our chosen halo is the second most mas- 
sive cluster in the ACDM simulation of Kauffmann et al 
(1999). We analyse its structure in the original simulation 
and in two higher resolution rcsimulations. In the collision- 
less case these are the lowest resolution members of a set 
of four resimulations carried out by Springel et al (2000a) 
using similar techniques to those of NFW. Details may 
be found there and in Springel ct al(2000b). These colli- 
sionless resimulations use GADGET as an N-body solver, 
whereas our collisional resimulations start from identical 
initial conditions but use the code's Smoothed Particle Hy- 
drodynamics (SPH) capability to solve the fluid equations. 
The SPH method regards each simulation particle as a 
"cloud" of fluid with a certain kernel shape. These clouds 
interact with each other over a length scale which is de- 
termined by the local density and so varies both in space 
and time. The basic parameters of our simulations are 
tabulated in Table 1, where Nt t 1S the total number of 
particles in the simulation, N n jg n the number of particles 
in the central high-resolution region, m p is the mass of each 
high-resolution particle, and l s stands for the gravitational 
softening length. Our cosmological model is flat with mat- 
ter density £l m = 0.3, cosmological constant 17a = 0.7 and 
expansion rate H = 70km~ 1 Mpc~ 1 . It has a CDM power 
spectrum normalised so that as = 0.9. The virial mass of 
the final cluster is M 2 oo = 7.4 x 10 14 /i _1 M Q , determined 
as the mass within the radius i?200 = 1.46/i _1 Mpc where 
the enclosed mean 



Table 1 
Simulation parameters 



Run 


iVtot 




m p (h- L M e ) 


l s (h I kpc) 


SO 


3.2x10° 


0.2x10° 


1.4 xl0 iu 


30 


SI 


3.5x10° 


0.5x10° 


0.68 xlO 1 " 


20 


S2 


5.1x10° 


2.0x10° 


0.14 xlO IU 


3.0 



3. RESULTS 

On scales larger than the final cluster, the matter dis- 
tribution in all our simulations looks similar. This is no 
surprise. The initial conditions in each pair of simulations 
are identical, so particle motions only begin to differ once 
pressure forces become important. Furthermore the initial 
perturbation fields in simulations of differing resolution are 
identical on all scales resolved in both models, and even 
SO resolves structure down to scales well below that of 
the cluster. As is seen clearly in Figure 1, a major dif- 
ference between the collisional and collisionless models is 
that the final cluster is nearly spherical in the former case 
and quite elongated in the latter. The axial ratios deter- 
mined from the inertia tensors of the matter at densities 
exceeding 100 times the critical value are 1.00:0.96:0.84 
and 1.00:0.72:0.63 respectively. Again this is no surprise. 
A slowly rotating fluid body in hydrostatic equilibrium 
is required to be nearly spherical, but no such constraint 
applies in the collisionless case (see also Miralda-Escude 
2000). 




Fig. 1. — The projected mass distribution in our two highest res- 
olution simulations. The collisionless case (S2) is on the top and 
the fluid case (S2F) is on the bottom. The region shown is a cube 
of 15/i -1 Mpc on a side. 



In Figure 2 we show circular velocity profiles for our 
simulations. These are defined as V c (r) — ^/GM(r)/r, 
where M(r) is the mass within a sphere radius r; they are 
plotted at radii between 21 s and 5i?200- They agree rea- 
sonably well along each sequence of increasing resolution, 
showing that our results have converged numerically on 
these scales. Along the fluid sequence the profiles resem- 
ble the collisionless case over the bulk of the cluster. In the 
core, however, there is a substantial and significant differ- 
ence; the fluid cluster has a substantially steeper central 
cusp. The difference extends out to radii of about 0.5i?2oo 
and has the wrong sign to improve the fit of CDM ha- 
los to published rotation curves for dwarf and low surface 
brightness galaxies. 
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Fig. 2. — Circular velocity profiles for our cluster simulations, 
each normalized to its own i?200 an d V^oo- These are plotted be- 
tween twice the gravitational softening and 5i?200- The collisionless 
sequence is plotted using dashed lines and the fluid sequence using 
solid lines. 

(Note that in the fluid case we expect small halos to 
approximate scaled down but slightly more concentrated 
versions of cluster halos, as in the collisionless case studied 
by Moore et al (1999a); this scaling will fail for interme- 
diate cross-sections because the ratio of the typical mean 
free path to the size of the halo will increase with halo 
mass.) 

In Figure 3 we compare the level of substructure within 
-R200 hi our various simulations. Subhalos are identified 
using the algorithm SUBFIND by Springel (1999) which 
defines them as maximal, simply connected, gravitation- 
ally self-bound sets of particles which are at higher local 
density than all surrounding cluster material. (Our SPH 
scheme defines a local density in the neighbourhood of ev- 
ery particle.) Using this procedure we find that 1.0%, 
3.4% and 6.7% of the mass within i? 2 oo is included in 
subhalos in SO, SI and S2 respectively. Along the fluid 
sequence the corresponding numbers are 3.0%, 6.4% and 
3.1%. The difference in the total amount results primarily 
from the chance inclusion or exclusion of infalling massive 
halos near the boundary at i?2oo- In Figure 3 we show 
the mass distributions of these subhalos. We plot each 
simulation to a mass limit of 40 particles, corresponding 
approximately to the smallest structures we expect to be 
adequately resolved in our SPH simulations. Along each 
resolution sequence the agreement is quite good, showing 
this limit to be conservative. For small subhalo masses 
there is clearly less substructure in the fluid case, but the 
difference is more modest than might have been antici- 
pated. 

4. SUMMARY AND DISCUSSION 

An interesting question arising from our results is why 
our fluid clusters have more concentrated cores than their 
collisionless counterparts. The density profile of an equi- 
librium gas sphere can be thought of as being determined 
by its Lagrangian specific entropy profile, i.e. by the func- 



tion m(s) defined to be the mass of gas with specific en- 
tropy less than s. The larger the mass at low specific 
entropy, the more concentrated the resulting profile. Thus 
our fluid clusters have more low entropy gas than if their 
profiles were similar to those of the collisionless clusters. 
The entropy of the gas is produced by a variety of accretion 
and merger shocks during the build-up of the cluster, so 
the strong central concentration reflects a relatively large 
amount of weakly shocked gas. 
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Fig. 3. — The total number of subhalos within -R200 is plotted as 
a function of their mass in units of M200. Dashed and solid lines 
correspond to the collisionless and fluid cases respectively. Results 
for each simulation are plotted only for halos containing more than 
40 particles. 

We study gas shocking in our models by carrying out 
one further simulation. We take the initial conditions of 
SI and replace each particle by two superposed particles, 
a collisionless dark matter particle containing 95% of the 
original mass and a gas particle containing 5%. These two 
then move together until SPH pressure forces are strong 
enough to separate them. The situation is similar to the 
standard 2-component model for galaxy clusters except 
that our chosen gas fraction is significantly smaller than 
observed values. 

In this mixed simulation the evolution of the collision- 
less matter (and its final density profile) is almost identical 
to that in the original SI. This is, of course, a consequence 
of the small gas fraction we have assumed. In agreement 
with the simulations in Frenk et al (1999) we find that 
the gas density profile parallels that of the dark matter 
over most of the cluster but is significantly shallower in 
the inner ~ 200/i~ 1 kpc. Comparing this new simulation 
(SIM) with its fluid counterpart (S1F) we find that in 
both cases the gas which ends up near the cluster centre 
lay at the centre of the most massive cluster progenitors 
at z = 1 ~ 3. In addition it is distributed in a similar 
way among the progenitors in the two cases. In Figure 4 
we compare the specific entropy profiles of the cluster gas. 
These are scaled so that they would be identical if each 
gas particle had the same shock history in the two simu- 
lations. Over most of the cluster there is indeed a close 
correspondence, but near the centre the gas in the mixed 



4 



simulation has higher entropy. (This corresponds roughly 
tor < lOO/i^kpc.) 
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Fig. 4. — We plot Lagrangian specific entropy profiles for the gas 
fluid simulation (S1F: crosses) and for the mixed simulation (SIM: 
open circles). In each case m(s) is given in units of the individual 
gas particle mass, m g , and the specific entropy of a particle is de- 
fined as ln(m 9 T g 1 - 5 /p 9 ). The arrows indicate where the timescale 
t2b for 2-body heating of the gas by encounters with dark matter 
particles (see equation (5) of Steinmetz & White (1997)) is 0.1, 1, 
and 10 times the age of the Universe. For each s we calculate t 2 b at 
the radius where the median specific entropy equals s. The dashed 
line with open squares is an "entropy" profile for SI calculated by 
using the SPH kernel to calculate the density and velocity disper- 
sion in the neighborhood of each particle, and then converting from 
velocity dispersion to temperature using the standard relation for a 
perfect monatomic gas. 

As Figure 4 shows, this is partly a numerical artifact; 
the two entropies differ only at radii where two-body heat- 
ing of the gas by the dark matter particles is predicted to 
be important in the mixed case. (The effect is absent in 



the pure fluid simulation.) The weaker shocking in the 
fluid case is evident from the equivalent "entropy" profile 
of SI in Figure 4. This lies between those of the two fluid 
simulations, and in particular significantly above that of 
S1F in the central regions. 

In conclusion the effective heating of gas by shocks in 
the fluid case is similar to but slightly weaker than that in 
the mixed case. This is presumably a reflection of the fact 
that the detailed morphology of the evolution also corre- 
sponds closely. The difference in final density profile is a 
consequence of three effects. In the mixed case the gas is 
in equilibrium within the external potential generated by 
the dark matter, whereas in the pure fluid case it must 
find a self-consistent equilibrium. In addition the core gas 
is heated by two-body effects in the mixed case. Finally in 
the pure fluid case the core gas experiences weaker shocks. 

Overall our results show that in the large cross-section 
limit collisional dark matter is not a promising candidate 
for improving the agreement between the predicted struc- 
ture of CDM halos and published data on galaxies and 
galaxy clusters. The increased concentration at halo cen- 
tre will worsen the apparent conflict with dwarf galaxy 
rotation curves. Furthermore, clusters are predicted to be 
nearly spherical and galaxy halos to have similar mass in 
substructure to the collisionless case, although with fewer 
low mass subhalos. Intermediate cross-sections would lead 
to collisional behaviour in dense regions and collisionless 
behaviour in low density regions with a consequent break- 
ing of the approximate scaling between high and low mass 
halos. The resulting structure may not lie between the two 
extremes we have simulated. Self-interacting dark matter 
might then help resolve the problems with halo structure 
in CDM models, if indeed these problems turn out to be 
real rather than apparent. 

SW thanks Jerry Ostriker and Mike Turner for stimu- 
lating discussions which started him thinking about this 
project. 
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